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j^fy Abstract 

With the goal of measuring localization in disordered interacting systems, we examine the finite-size scaling of the geometrically- 
averaged density of states calculated from the local Green's function with finite energy resolution. Our results show that, unlike in 
| I a simple energy binning procedure, there is no limit in which the finite energy resolution is irrelevant. 
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Significant disorder, for example due to non-stoichiometric 
doping, is a feature of many materials, and may in fact 
play a key role in many interesting behaviors, especially 
when combined with electronic correlations. The study of 
systems with both disorder and interactions has a long 
history, pQ but it has been difficult to access the regime in 
which both are strong, where the physics of both Anderson 
localization and the Mott gap are significant. The tech- 
nique of dynamical mean field theory 2J (DMFT) has been 
very successful in the study of systems with strong inter- 
actions, and there has recently been interest in applying 
it to disordered systems. [3|4j DMFT calculates the local 
Green's function, G(r, r, to + i"f). How can information on 
localization be drawn from this? 

The geometrically-averaged density of states (GADOS) 
appears to be ideally suited to the task. This is the geomet- 
ric average over sites of the local density of states (LDOS) 
which is obtained directly from the local Green's function: 
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which is exact in the limit 7 — > + . In an infinite system, 
p g {to) at the Fermi level is an order parameter for the An- 
derson transition, being nonzero when states at the Fermi 
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level are extended and zero when they become localized. 
This is because when wavefunctions at a given energy be- 
come localized in real space, the LDOS at that energy be- 
comes discrete: localized wavefunctions at the same energy 
are well separated in position so the LDOS at that energy 
goes to zero far from the wavefunction centers. 

To consider specific disorder configurations, it is neces- 
sary to work with finite size systems, and in this case finite- 
size scaling analysis of the GADOS must be performed. 
This procedure is well studied in non-interacting systems 
for which single-particle wavefunctions can be calculated, 
providing perfect energy resolution. In a numerical calcu- 
lation on an interacting system, however, the Green's func- 
tion is obtained for finite 7. We have explored the impact 
this has on the scaling of the GADOS in the context of a 
noninteracting model 
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with site energies, ej, distributed with probability -P(ei) = 
6(W/2 - |ei|)/W. The width of the potential distribution, 
W, parameterizes the strength of the disorder. The model 
is solved exactly by numerical diagonalization using the 
LAPACK libraries. In earlier work, [5] we calculated the 
local density of states within a finite energy resolution 7 
using a simple binning procedure: 
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Each state which falls in the energy window is assigned 
equal weight. In this case, we demonstrated that the GA- 
DOS scales as a function of an effective length scale: L e t f — 
1/(1/1/ + 1 /£ 7 ) where £ 7 acts as an effective inelastic mean 
free path, limiting the relevant system size. t 1 can also be 
seen as the system size divided by the average number of 
states in a bin: n ~ 7/ A where in one dimension A = 
l/(pL) is the average energy spacing between states and 
p{uo) is the density of states. 

Here we model more accurately the way the Green's func- 
tion introduces finite energy resolution: 
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All results shown are at the band center on one-dimensional 
systems with disorder strength W = 1, for which all states 
are localized. 
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Fig. 1. The geometrically-averaged density of states at the band 
center in a one-dimensional system with disorder strength W = 1. All 
points represent averages over at least 1000 disorder configurations. 
L is the system size, and £ e // = L/(pA) where p is the density 
of states, (a) Square energy bin (Eqn. (|4jl) with constant width 7. 
(b) Lorentzian energy bin (Eqn. JSJ) with constant 7. (c) and (d) 
Lorentzian energy bin with 7 = A/L. 

Fig. 1(a) shows the GADOS as a function of system size 
for three non-zero values of 7 using the square binning in 
Eqn. (g]) reported earlier. [5j Fig. 1(b) shows the GADOS 
as a function of system size for the same 7 values using 
the Lorentzian binning in Eqn. (0. Two features are im- 
mediately apparent: the lack of exponential decay and the 
overall inflation of the GADOS. With the Lorentzian bin, 
all states contribute to the sum in Eqn. ([5]) with rapidly 
decreasing weight away from to. Despite the small weight 
assigned to them, the many states outside the peak of the 
Lorentzian have a dramatic impact on the GADOS. This 
is because when these states are centered, in real space, at 
positions with otherwise small LDOS they can significantly 
increase the LDOS. Even a modest increase in the mini- 
mum values of the LDOS causes a dramatic increase in the 
GADOS. 



A second important distinction between the two binning 
schemes arises in the limit 7 < A. In this limit, for each 
u> there will be disorder configurations for which (i) in the 
square bin case there are no states in a bin and (ii) for 
the Lorentzian bin case no state has significant weight. In 
the case of the square bin, we removed these bins from 
the average. However, there is no systematic way to do 
this in the Lorentzian case. The inclusion of these states 
dramatically lowers the configuration averaged LDOS and 
hence the GADOS even below the 7 = value. 

It appears to be possible to account for the behavior of 
p g with Lorentzian binning within the framework of an ef- 
fective length; however, a different form of L e // is required. 
In the square bin case, in the limit of small 7, L e ff = L, 
because by construction there is always at least one state 
in each bin. In the Lorentzian bin case, the average number 
of states in a bin may be less than one. £ e // is now simply 
equal to L/fi = .LA/7 = 1/(^7). This is consistent with 
the nearly size independent results when a fixed value of 
gamma is used. Furthermore, Fig. 1(c) shows the GADOS 
as a function of system size using bin widths 7 = A/L. Fig. 
1(d) presents the same data plotted versus 
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showing scaling. The key point is that there is no limit in 
which the finite energy resolution is irrelevant. 

In conclusion, there appears to be no finite-size scaling 
of the GADOS calculated from the local Green's function 
with fixed energy resolution, and the scaling which occurs 
when 7 cx 1/L arises at least in part from physics unrelated 
to localization. We suggest that other measures of localiza- 
tion, such as the inverse participation ratio, which do not 
share the extreme sensitivity to small values of the LDOS 
may be more effective than the GADOS in interacting sys- 
tems, although changes to the scaling due to finite energy 
resolution will still need to be addressed. 
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